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We study theoreticahy the photoabsorption spectra of an attosecond XUV pulse by a laser- 
dressed atomic system. A weak XUV excites an autoionizing state which is strongly coupled to 
another autoionizing state by a laser. The theory was applied to explain two recent experiments 
[Loh, Greene, and Leone, Chem. Phys. 350, 7 (2008); Wang, Chini, Chen, Zhang, Cheng, He, 
Cheng, Wu, Thumm, and Chang, Phys. Rev. Lett. 105 143002 (2010)] where the absorption 
spectra of the XUV lights were measured against the time delay between the laser and the XUV. 
In another example, we study an attosecond pulse exciting the 2s2p(^P) resonance of helium which 
is resonantly coupled to the 2s^(^ S) resonance by a moderately intense 540 nm laser. The relation 
between the photoabsorption spectra and the photoelectron spectra and the modification of the 
transmitted lights in such an experiment are analyzed. The role of Rabi flopping between the 
two autoionizing states within their lifetimes is investigated with respect to the laser intensity and 
detuning. 

PACS numbers: 32.80.Qk,32.80.Zb,32.80.Fb,42.50.Gy 



I. INTRODUCTION 



Autoionization in an atomic system is understood as 
the decay of a bound state to one or more degenerate 
continua. It has long been observed in photoabsorption, 
photoion, or photoelectron spectra where resonances are 
characterized by their detailed shapes as given by Fano's 
configuration interaction theory In such energy- 

domain measurements the lifetimes are deduced from the 
energy widths of the resonances. With the development 
of ultrafast technologies in the past decade, in particu- 
lar, the emergence of single attosecondpulscs (SAPs) in 
the extreme ultraviolet (XUV) region it has become 
possible to ask whether one can use a very short pulse 
to populate an autoionizing state and watch its decay 
within the lifetime. This issue has been addressed the- 
oretically In Ref. we proposed a pump-probe 
scheme using two attosecond XUV pulses where the time 
evolution of the energy profile of the autoionizing wave 
packet initiated by the pump is followed by another at- 
tosecond pulse within the decay lifetime. However, such 
measurements are not yet possible because the available 
attosecond pulses are still too weak. Instead, in actual 
experiments today, attosecond pulses are used only to 
excite an autoionizing state which is then "probed" by 
an intense infrared (IR) field 0,[1|- However, the XUV 
and the IR pulses often overlap in time, and thus it is 
not appropriate to consider the IR only as a probe pulse, 
especially when the IR is rather intense. Instead, it is 
preferable to consider such XUV-I-IR experiments as pho- 
toabsorption or photoionization of a laser-dressed atom. 
If photoionization occurs in the energy region where the 
electron continuum spectrum is featureless, i.e., with- 
out bound states or resonances, then the presence of the 
IR can be incorporated using the approximate "streak- 
ing" model 0. A number of such experiments using 



SAPs and attosecond pulse trains (APTs) have been per- 
formed [l0l - [l3| and explained based on the attosecond 
streaking model. However, bound states and autoioniz- 
ing states are prevalent in atomic and molecular systems. 
The presence of an intense laser is likely to resonantly or 
near-resonantly couple XUV-initiated states with other 
bound or autoionizing states. In the limit of long pulses, 
such strongly coupled systems have been widely inves- 
tigated in the past few decades as a model three-level 
atomic system, and interesting phenomena such as elec- 
tromagnetically induced transparency (EIT ) [li.fisl . co- 
herent population trapping [3] and others jl5l 17l| have 
been widely investigated. Moreover, they are generally 
perceived as an all-optical control of several characteris- 
tics of the light beam propagating in the medium, such as 
the speed, absorption, storage, and retrieval of light [Tsl - 
[20j . In most of these studies, the coupling laser and the 
probe laser (or the pump laser) are limited to visible and 
infrared lights. However, the weak probe (or pump) laser 
can be extended to the high-energy spectral region from 
XUV to X-rays. These light sources excite inner-shell 
electrons or multiply excited states of an atom while the 
laser strongly couple these states to other nearby atomic 
states. In fact, EIT-like phenomena have been reported 
in such systems [2l| and EIT-like theory has also been 
extended to the systems where the two strongly coupled 
states are autoionizing states [22| - [25j . For this case with 
long pulses, one can measure both absorption and elec- 
tron spectra by scanning the wavelength of the probe 
pulse. 

With the emergence of SAPs, we recently general- 
ize the three-level autoionizing systems for the situation 
where both the XUV and the laser are short pulses [1^ . 
Specifically, we modeled the experiment by Gilbertson et 
al. where a 140-attosecond pulse with central energy 
of about 60 eV is used to excite the 2s2p{^P) doubly 
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excited state of helium in the presence of a moderately 
strong few-femtosecond IR laser. The IR laser can couple 
2s2p to the 2p'^(^S) doubly excited state. This appears to 
be a typical three-level system where one can investigate 
the autoionization dynamics of 2s2p within its lifetime 
(17 fs), or more precisely, the use of IR to control autoion- 
ization dynamics by changing the time delay between the 
two pulses. Unfortunately, due to the large IR ionization 
rate from the higher 2p^ state, the strong IR coupling in 
this case mainly provides an efficient pathway to remove 
electrons from 2s2p, thus depleting the autoionized elec- 
trons which would be observed in the time delay spectra. 
In order to observe the change of the resonance profile, 
spectral resolution much better than 0.1 eV is needed. 
In the measurement the electron spectrometer has a 
resolution of 0.7 eV. Thus, the lifetime information has 
been retrieved, but no information on the autoionization 
dynamics has emerged from this experiment. 

To achieve better spectral resolution, measuring pho- 
toabsorption is more favorable. An XUV-I-IR experiment 
was reported by Wang et al. [g^l where the absorption 
spectra of the laser-coupled autoionizing states in argon 
were measured against the time delay. In this paper, we 
extend our earlier theory p6j to obtain absorption spec- 
tra versus time delay to explain the results in Ref. [27|. 
Practically, the effect of the light propagation in the gas 
medium has to be incorporated [28|, which we will ad- 
dress in the future. The present theory is also applied 
to the transient absorption spectra with pulse trains re- 
ported by Loh et al. [21|, where the lengths of the two 
light pulses, 30 and 42 fs, were longer than the lifetimes 
of the autoionizing states of interest. To fully utilize the 
laser coupling, we suggest a 540 nm laser to couple reso- 
nantly 2s2p and 2s^{}S) in helium. In such a setup, these 
two autoionizing states are very weakly ionized by the 
laser and can possibly undergo one or more Rabi oscil- 
lations within their lifetimes. By changing the intensity 
and the wavelength of the laser, there are rich possibili- 
ties to manipulate autoionization, and the effects can be 
measured in either the absorption or the electron spec- 
tra. Unlike long pulses, the large bandwidth of an at- 
tosecond pulse covers the whole autoionization spectrum 
in a single measurement for each fixed time delay. Such 
manipulation would epitomize true attosecond control of 
electron dynamics. In the meanwhile, the modification 
of the XUV absorption, if incorporated with appropri- 
ate macroscopic conditions of the medium, would sub- 
stantially alter the transmitted pulse going through the 
medium. This in turn would offer a tool to shape ultra- 
short light pulses, which serves in the theme of optical 
coherent control j29l - l32| . 

In Sec. ini our model for the total wavefunction is in- 
troduced. It is generalized from the earlier models for 
long pulses [l^ to short pulses. Then we look at ab- 
sorption in terms of the frequency-dependent response 
function, which is defined as the energy-differential prob- 
ability density of the absorption through dipole interac- 
tion [28|. In Sec, mil the calculated spectra are compared 



with the available experiments (2lJ,|27|, and in a test case, 
the absorption and the electron spectra are compared and 
discussed. The conclusions are drawn in Sec. IIVI Atomic 
units (c = ?i = e = 1) are used throughout Sec. [Tll and 
electron Volts (eV) for energy and attoseconds (as) or 
femtoseconds (fs) for time are used in the rest of this 
paper unless otherwise specified. 



II. THEORY 
A. Total wavefunction 

The total wavefunction of a laser-assisted three-level 
autoionizing system has been formulated in our previous 
report [2^. Consider a model atom consisting of the 
ground state \g) and two autoionizing states \a) and 
which are degenerate with the continuum states \Ei) and 
I £'2) respectively. The atom is exposed to an XUV field 
that couples \g) to the first resonance and an IR field that 
couples the two resonances. The pulses are in arbitrary 
envelopes and linearly polarized in the same direction. 
The total Hamiltonian is H{t) = Ha + Hi{t) where Ha 
is the atomic Hamiltonian and Hj{t) is the interaction 
between the atom and the light. The field is treated as 
an external disturbance to the system where the field's 
own Hamiltonian is excluded. The total wavefunction 
can be written in the configuration space as 



|vl/(i))=e--«'c,(t)|5) 



da{t)\a)+ J dEAt)\Ei)dEi 



dbit)\b) + / dE,it)\E2)dE2 



(1) 



where Eg is the ground state energy, and Ex = i-^x + Eg 
and El = ujl + tJx + Eg are the atomic energies popu- 
lated by the pulses. In principle, the energy integrations 
are over infinite range, but practically only the parts near 
the resonances arc infiucnccd by the strong coupling, thus 
worth our concern. Since the fast oscillating parts are 
factored out in the exponential terms in Eq. ([1]), the co- 
efficients Cg(i), dait), dEi(t), db{t), and dE^it) are all 
slowly varying functions of time. This form is valid when 
the field is weak to moderate such that the dipole tran- 
sitions between \g) and \a) and between |a)and \b) are 
the major dynamics being considered. This means mul- 
tiphoton effect or tunneling ionization is to be treated 
simply as additional depletions of the bound states \a) 
and 1 6), where the rates can be estimated by the theory 
developed by Ammosov, Delone, and Krainov (ADK the- 
ory) [33| or the theory developed by Pcrclomov, Popov, 
and Terent'ev (PPT theory) [sl]. 

Substituting the expansion of |^'(i)) into the time- 
dependent Schrodingcr equation, first-order coupled 
equations for the five coefficients are found, where two 
coefficients are continuous. As described in 1261 . we first 
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apply adiabatic elimination to dE^it) and dsiit) to ob- 
tain the "first iteration" of these two functions of en- 
ergy and time. They are then plugged into the other 
three coupled equations for Cg{t), da{t) and db{t) (See 
Eqs. (12-14) in Ref. [1^), where these bound-state coef- 
ficients can be solved numerically as functions of time. 
Next, in the "second iteration", correcting the adiabatic 
elimination, the continuum-state coefficients arc solved 
with the original coupled equations (See Eqs. (18) and 
(19) in Ref. Thus, with ah the coefficients in Eq. ([T]) 

calculated, |vE'(t)) is obtained. In our model, the dipole 
matrix elements and the Fano parameters are assumed to 
be independent of electron energy, and they can be ob- 
tained from energy-domain experiments or from separate 
structure calculations. 

The total wavefunction |^'(t)) can also be expressed in 
the eigenstate basis as 

\^{t)) = e-^«*c,(t)|5) +6-^-*! 4")(0|V'i"^)di? 

/ c^^}m^^})dE\ (2) 



where the eigenstates IV'^'') and {ip^h composed of 
I a) and \Ei) and of \b) and [£'2) respectively, as described 
in Fano's theory [l|. When the field is over, and 
\ipE^) are stationary states of the atomic Hamiltonian, 

and the coefficients c'"^\t) and c-^ {t) are constant of 
time other than a factor of oscillating phase. The pho- 
toelectron spectra of the two resonances are calculated 
as P^'^^E) = \c^E\tf)? and P^^\E') = |c|^(t/)P, re- 
spectively, where i/ is taken after the pulses are over. In 
this paper, the time delay is defined between the pulse 
peaks, and it is positive if the IR comes after the XUV. 



B. Absorption 

The frequency-dependent response function S{uj) for 
photoabsorption of an atom exposed to a finite-duration 
external field has been described recently by Gaarde et 
al. [2^. Considering an atom in an oscillating electric 
field, the response function S{La) is defined by 



AJ7= / ujS(uj)duj 
Jo 



(3) 



where A[/ is the total energy absorbed by the atom from 
the field. The response function represents the probabil- 
ity density for the light of frequency uj to be absorbed. 
The value of S{u!) is positive for absorption and nega- 
tive for emission, for a positive to. In our model, the 
light-atom interaction term in the Hamiltonian is given 
by —D{t)E{t) where D{t) is the electric dipole, and the 



total absorbed energy through the dipole interaction is 
dD{t) 



AU = 



dt 



-E{t)dt 



(4) 



= -2/ ujlm[D{uj)E*{u)]dw, (5) 
Jo 

where D{uj) and E{ijj) are Fourier transforms of D{t) and 
E{t) respectively. The convention of Fourier transform 
for an arbitrary finite function /(i), from the time do- 
main to the frequency domain, is 



1 



'f{t)dt. 



(6) 



Note that in Eq. ([5]) we have implemented the fact that 
D{t) and E{t) are all real functions so that D*{lu) = 
b{—oj) and E*{ijj) = E{—oj). The final form of the re- 
sponse function, as extracted from Eq. ^ and Eq. 



S{u) = -2lu\[iD{uj)E* {uj)]. 



(7) 



As mentioned earlier, 5(0;) is the absorption probability 
density, which has the same dimension of the photoelec- 
tron profile P{E) defined in Sec. Ill Al The absorption 
cross section a{uj) is related to ^(a;) by 



a{uj) 



\EH\ 



2 ' 



(8) 



where a is the fine structure constant. 

With the total wavefunction in Eq. ([J), the dipole mo- 
ment D{t) is given by 

D{t) ^ {^{tM^{t)) (9) 
= [e^"^*ux(t) + e*"^*Ui(t) + c.c] (10) 

where ^ is the electric dipole operator (/x = —ez), and 
ux(t) and UL{t) are 



ux{t) = cg{t) j M^^K^^^*{t)dE (11) 
UL{t) ^ I c^;^\t) f Mi%c^^}*{t)dE'dE, (12) 



where Af^?) ^ (^^Vlff) and mI,% ^ (^i'VlV'i"^) 
are the dipole matrix elements between the eigenstates. 
With Fano's theory [H, M^"^ and M^j^^ 

can be decom- 
posed to the dipole matrix elements in the configuration 
space, which are taken as input parameters in our model. 
Following Eq. (|10p . the D{uj) function is 



^/2^ J- 



s'("^-")*ux(i) + e''("^-'^)*Wi(i)]dt 

(13) 

where the rotating wave approximation has been applied, 
provided that Et, > Ea- U Ef, < Ea, u>l will be negative. 
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and e*'^^*UL(i) should be substituted by e '■'^^*u*^{t) in 
Eq. ([T^ . Since the contributions of the two Ught pulses 
to D{ui) are quite separate in frequency in the current 
setup, we can write the response function involving ux{t) 
and Exit) as S'x(w), and that involving UL{t) and EL{t) 
as Sl{^)- 

In the following calculations, the photoclectron profiles 
for the two resonances are P^"^ (E) and P^''-' (E) defined in 
Sec. Ill Al and the photoabsorption profiles for XUV and 
IR arc S'x(w) and Sl{^), respectively. If the density of 
the atomic gas is sufficiently low, the profiles given above 
will be linear to the actual measured electron and light 
signals, which can then be assessed at an independent 
stage. 



III. RESULTS AND DISCUSSIONS 

A. Photoabsorption spectra of Ar with single 
attosecond pulses 

We now apply the present theory to analyze the exper- 
iment by Wang et al. [l^] where an XUV SAP was used 
to excite argon atoms to the 3s3p^4p autoionizing state 
in the presence of an intense IR pulse which strongly cou- 
ples this state with the 3s3p^4d autoionizing state. The 
weak XUV pulse had duration (tx) of 140 as and pho- 
ton energy (uix) covering 20 to 40 cV. The IR pulse had 
wavelength (Al) of 750 nm [lol = 1-65 eV) and dura- 
tion (tl) of 7 fs. Two IR peak intensities (/l), 0.5 and 
1.0 TW/cm^ were independently applied. The experi- 
ment measured light transmission of the XUV for differ- 
ent time delays. The experimental results for the two IR 
intensities are shown in Fig.[TJa-b). The main conclusion 
of the experiment is that the transmission peak at the 
Fano resonance is depressed when the two pulses overlap, 
and when the IR lags behind the XUV, the main peak re- 
covers gradually as the two pulses moved away from each 
other on the scale of the decay lifetime of 3s3p^4p. This 
result serves to reveal autoionization in the time domain. 
Note that the delay decreases along x-axis in Fig. [U so 
the positive delay (XUV first) arc on the left of the plots, 
which follows the graphing convention in Ref. [27| . 

Considering a dilute target gas, the transmission spec- 
trum can be written as 



'> 27-0 



T(u!) = Tq{uj) cxp[— p(T(a;)L] 



(14) 



where T(uj) and To(aj) are the light intensity going into 
and coming out of the medium respectively, p is the gas 
density, o'Ilu) is the absorption cross section for a single 
atom, and L is the length of the light path through the 
medium. The approximation « 1 — a; can be fur- 
ther applied for a low-density gas, where the final mea- 
surement is linearly related to the single-atom response. 
In this paper, the calculations are made only for single 
atoms, where the information of the medium, which is 
often unavailable, is not needed. 



26.0 26.0 

40 30 20 10 -10 40 30 20 10 -10 



? 27.0 





26.0 25.0 

40 30 20 10 -10 40 30 20 10 -10 



(f) 


s 







1.0 

^ 0.8 

0.6 
H 0.4 
0.2 
0.0 



40 30 20 10 10 
Delay [fs] 



40 30 20 10 -10 
Delay [fs] 



FIG. 1: (Coler online) Time-delayed transmission spectra of 
the 3s3p^4d resonance in Ar for two IR intensities. The signals 
are normalized to 1. The measurements shown in (a-b) are 
taken from Ref . 12711 . The simulations in (c-d) are based on the 
model in Ref. [23] and updated by Zhang [Sgl. The present 
calculations are shown in (e-f). All left panels are for II = 
0.5 TW/cm^ and all right panels are for II = 1.0 TW/cm^. 



In Ref. [27|, the simulation was carried out and re- 
cently updated by Zhang [33. His new results arc shown 
in Fig. (Ijc-d). Unlike in Ref. [23], the spectrum is now de- 
fined as the response function as in Eq. ([7]). The dipole 
moment therein is calculated with the total wavefunction 
of the earlier model in Refs. [23l425j . which is equivalent 
to our model with the exact bound-state coefficients and 
the first iteration of the continuum-state coefficients (See 
Sec. Ill A} . As seen in Fig.jlja-d), the main discrepancies 
between the measurements and the updated calculations 
by Zhang are as follows. First, near to = 10 fs where the 
two pulses overlap, the calculation in (d) shows strong 
peak splitting, but the measurement in (b) has only the 
upward-curving branch. The same discrepancy is also 
seen between (a) and (c). Second, in the measurements, 
the resonance peak is depleted much more strongly at 
the high intensity than at the low intensity, while the 
simulations in (c-d) do not have this clear contrast. 

The results of the present calculation are shown in 
Fig. dfe-f). We include the tunneling ionization rates 
for the autoionizing states estimated by the ADK model 
with the parameters fit to the experiment. For 3s3p®4p 
and 3s3p^4(i, C; = 1.6 and 0.7, and a = 4 and 5, respec- 
tively, where C; is the ADK coefficient in Ref. [3j] and a 
is the correction parameter in Ref. [3^ . Our calculations 
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differ from tfie previous ones, as shown in Fig.[TJc-d) ver- 
sus (e-f), in the following ways. First, the resonance peak 
is heavily smeared by the IR ionization near the overlap 
of the pulses. This is especially obvious in Fig. [TJf). Sec- 
ond, the upper branch of the split peaks tilts upward 
more than that of the previous simulation, and is closer 
to what the measurement suggests. Third, below the res- 
onance energy, there are negative fringes in the spectra 
curving downward with the delay, which are not seen in 
the previous report. By comparing the two models with 
the experiment, it is obvious that the most important 
improvement by the present one is the inclusion of the 
IR ionization. The one more iteration of the continuum- 
state coefficients, which is newly developed in our model, 
is not essential in this case. Some other detailed features 
of the result, such as the fringes, are unexplained so far. 
On the other hand, they may be washed out in measure- 
ments where in principle, signals should be integrated 
over a range of intensities. 
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FIG. 2: The transient absorption spectra (upper panels) 
and the photoabsorption spectra (lower panels) reported in 
Ref. [SJ] (left panels) and calculated by the present model 
(right panels) for various time delays. The experimental data 
reported in Ref. [2j| is shown in (a) in dots. 



B. Transient absorption spectra of He with 
attosecond pulse trains 

Electromagnetically induced transparency (FIT) [T^ . 
[Tsj has been traditionally studied as a time-independent 
phenomenon where a resonance state shows a depleted 
profile in the absorption spectra when an ac field strongly 
couples it to another excited state. While it has been 
carried out by long pulses and low-lying excited states, 
a recent work by Loh et al. [2]| consisting of measure- 
ments and simulations investigated the FIT effect ini- 
tiated by femtosecond pulses. It employed an 800 nm, 
42 fs, 14 TW/cm^ IR pulse to couple the 2s2p{}P) and 
2p^(^S') doubly excited states in helium, and a 30 fs, 
weak XUV produced by high-order harmonic generation 
to probe the IR-coupled system at different time delays. 
The simulation therein employed the model developed 
by Lambropoulos |22| configured for longer pulses. Our 
analysis for this experiment illustrates that the present 
model can be applied to broad-band attosecond pulses as 
well as longer femtosecond pulses. 

Figure [H^a) shows that the theoretical model in 
Ref. |2l| can reproduce the major features of the transient 
absorption spectra observed in the experiment. Transient 
absorption spectrum is defined by absorption spectrum 
subtracted by static spectrum. Since the IR is rather in- 
tense, in Ref. [2l| it was concluded that the resonance 
widths of 2s2p and 2p^ should be corrected by the IR 
ionization using the tunneling ionization rates specified 
in the paper. The transient spectra in Fig. [2ja) have 
also been convolved with a 0.18-eV Gaussian function to 
account for the spectrometer resolution. A post-pulse 
which has 33% of the peak IR intensity, at 88 fs after 
the main peak is taken into account. The post-pulse is 
responsible for the small peak at about 0.45 eV in the 
absorption spectra in Fig. [2l^b) and (d). 

In Fig. mb), the "raw" theoretical photoabsorption 



spectra from Ref. 121| are shown. Here, the "static" spec- 
trum means the XUV spectrum measured in the absence 
of the IR. All the profiles are calculated as the ionization 
rate versus XUV detuning in weak-probe limit, given by 
Eq. (1) in Ref. [2l|, which employs the formalism by Mad- 
sen et al. [l^l- In this formalism, the XUV pulse shape 
is not considered, and as shown in the spectra, the reso- 
nance peaks are not broadened by the XUV duration. In 
our model, we calculate the profiles as either the ground 
state depletion 1— |cg(t/)p or the total absorption proba- 
bility J S{uj)dLo (which are equivalent theoretically) ver- 
sus XUV detuning. Our calculation assumes the XUV 
peak intensity Ix = 10^° W/cm^, such that it is in the 
first-order perturbation regime. All the XUV and IR 
electric field parameters are included in our model. The 
60-meV bandwidth of the XUV pulse broadens the reso- 
nance peaks in Fig. [5Jd). 

The theoretical absorption spectra show that when the 
two pulses overlap, the 2s2p resonance peak splits into 
two peaks moving outward, corresponding to the Autler- 
Townes doublet; at the same time, the peaks are highly 
suppressed due to the IR ionization. A good agreement 
has been found between the experiment and both mod- 
els. The FIT effect by femtosecond pulses, as the goal of 
this study, can thus be controlled by overlapping the two 
pulses with different delays. However, because of the nar- 
row XUV bandwidth, the absorption is defined for each 
XUV photon energy, and it is not necessary to carry out 
the energy distribution in each measurement by means of 
the response function S{uj). In other words, our model is 
not advantageous over Lambropoulos' model here. How- 
ever, for the argon case in Sec. IIII Al the broadband at- 
tosecond XUV retrieves the whole absorption spectrum 
in just one measurement, which can be treated in the 
present model but not in Lambropoulos' model. 
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C. Electron emission and photoabsorption with 
single attosecond XUV pulses 

The outcome of an XUV+IR experiment is usually de- 
fined by the time-delayed electron or absorption spectra. 
In the case of weak XUV and strong IR, as in EIT, the 
interaction can be viewed as an XUV probing a strongly 
IR-coupled atomic system. For the autoionizing states 
considered here, the IR can also be viewed as a way 
to control or manipulate the autoionization dynamics. 
In the case of photoabsorption, the single-atom response 
studied in this work can be implemented in the light field 
traveling through the medium. It can then be regarded 
as a form of optical control or pulse shaping of the XUV 
by the IR |30l - [32 |. which will be our extended effort of 
this project in the near future. Both the electron and 
the absorption spectra are the observable consequence 
of the interplay between the electron dynamics and the 
optical fields. In this subsection we examine the relation 
between these two types of measurements. Before we pro- 
ceed, note that while the modification of the XUV by the 
IR has often been studied, the IR is also modified by the 
XUV and complimentary "shaping" should appear also 
in the IR spectrum. In the following, by changing the 
laser wavelength to 540 nm (it is now actually a visible 
light instead of an IR) to couple 2s2p{^P) and 2s'^(}S) in 
helium, we will study the absorption spectra of both the 
XUV and the laser, along with the electron spectra. The 
laser spectra shown therein are not only informative and 
interesting, but also assessed as realistically measurable. 

In this test case, an XUV pulse with duration tx = 
100 as, peak intensity Ix = 10^" W/cm^, and cen- 
tral photon energy tox = 60 eV is shone on helium 
to excite it to 2s2p at 60.15 eV as well as the "back- 
ground" lsep{^P) continuum states. Because of the wide 
bandwidth, other nearby doubly excited states, such as 
2s3p{^P) and 2s4p(^P), are also excited as well; they 
have much narrower widths and are weaker [37| . thus 
disregarded in the present study. In the meanwhile, a 
time-delayed laser pulse with wavelength = 540 nm, 
duration tl = 9 is, and peak intensity II = 0.7 TW/cm^ 
is applied on the system to strongly couples 2s2p to the 
2s'^{^S) doubly excited state at 57.85 eV which is below 
2s2p energetically. The single-ionization binding energy 
of 2s2p and 2s^ are 5.3 eV and 7.6 eV, respectively. With 
the laser applied, their tunneling ionization rates are neg- 
ligible, and we expect to see the effect of strong coupling 
between them. In the simulation, the dipole matrix ele- 
ments between Is^ and 2s2p and between 2s2p and 
are 0.038 and 2.17 a.u. respectively. The Fano parame- 
ters are F — 37 meV (lifetime is 17 fs) and q = —2.75 for 
2s2p, and F = 0.125 eV and q = 1000 for 2s^. The direct 
ionization path {lses\fi\2s2p), which involves transition 
of two electrons, is expected to be very weak compared 
to (2s^|/x|2s2j>), where the ratio is determined by the q- 
parameter. The calculated electron and absorption spec- 
tra are done with various time delays to between -15 and 
50 fs. 



The electron and absorption profiles are shown in 
Fig. [3l The electron emissions near the two resonances 
are shown in Fig. |3l(a-b). The photoabsorptions for the 
XUV and the laser near their central photon energies are 
shown in Fig. [3jc-d). Note that the energy scale for all 
four graphs are the same, which means the signal dis- 
tributions for these four calculations can be compared 
side-by-side. The electron emission of this test case has 
been reported in details previously by us [1^; the main 
task here is to investigate its relation to photoabsorption. 
Two observations can be drawn from Fig. [3l 




-10 10 20 30 40 50 -10 10 20 30 40 50 

Delay [Is] Delay [fs] 

FIG. 3: (Color online) Photoelectron spectrograms near the 
(a) 2s2p and (b) 2s^ resonances, and photoabsorption spec- 
trograms for the (c) XUV and (d) laser fields, for a fixed laser 
intensity and wavelength (See text). The negative signals 
plotted in blue and cyan colors in (c-d) represent photoemis- 
sion. 

First, the spectrogram (signal density versus energy 
and time delay) of the 2s2p resonance is almost the same 
in the electron and the absorption profile. However, the 
absorption probability is almost ubiquitously higher be- 
cause the electrons brought by the XUV photons from 
Is^ to 2s2p are redistributed between the 2s2p and 2s^ 
resonances via the laser coupling. In other words, all the 
electrons shown in Fig. EJa) and (b), coming out from 
2s2p and 2s^ respectively, have gone through the ab- 
sorption shown in Fig. |3l(c). This observation is further 
clarified by energy conservation of the system. Figure |4] 
shows the total energies absorbed by the XUV and the 
laser, as well as the total emitted photoelectron energies, 
as functions of time delay. At each time delay, the total 
energy absorbed by the XUV is calculated by integrating 
S{uj) over the 1-eV range centered at the resonance en- 
ergy 60.15 eV. Outside this region, the electrons are not 
involved in the laser coupling, thus not changing with 
the time delay. The total emitted electron energy of the 
2s2p resonance is calculated by integrating E\c^^\tf)\'^ 
over the same range. The total energy of laser absorp- 
tion is small due to its low photon energy and narrow 
bandwidth. The total emitted electron energy of the 2s^ 
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resonance is obtained by counting energy over the whole 
resonance width (Recall that q = 1000 indicates that the 
IR couphng does not generate the Ises continuum). Fig- 
ure S] shows that when the laser pulse appears just after 
the XUV, a considerable fraction of energy shifts from 
2s2p to 2s^. As the delay increases, this fraction drops 
because more electrons decay from 2s2p before the laser 
couples. The work of laser is to direct electrons to the 
decay paths of 2s2p or 2s^, which are separate by only 
2.3 eV. Energy conservation is obvious where the total 
absorbed energy from XUV and from laser equals the 
total emitted energy from 2s2p and from 2s^. 




-10 10 20 30 40 50 

Delay [Is] 



FIG. 4: The total absorbed and emitted energy across 1-eV 
spectral range over the resonances versus time delay. The 
solid curves are for electron emission and the dashed curves 
are for photoabsorption. 

The second observation is the strong correlation be- 
tween the electron profile of the 2s^ resonance and the 
absorption profile of the laser for the time delay between 
and 5 fs, as seen in Figs. [3Kb) and (d). The 2s^ profile 
shows a major peak for which the stimulated emission 
of laser drives the helium from 2s2p to 2s^ before 2s2p 
decays. The laser absorption profile has a negative peak 
of the same strength and in the same time-delay range 
as the positive peak of the 2s^ electron profile. In other 
words, the Rabi oscillation during and 5 fs is shown 
by the electron burst and the laser emission equivalently. 
However, as the delay increases, the absorption profile 
quickly becomes complicated. For each time delay, the 
interference fringes are shown in both the absorption and 
the emission parts of the laser spectra. No obvious ex- 
planations of such structures can be offered at this time. 



D. Intensity and wavelength dependence 

The spectrograms shown in Fig. |3] contain rich details 
of how the dressing laser controls the autoionization dy- 
namics and the absorption of the XUV light. Such details 
are observable only if electrons and photons are detected 
with high spectral resolution because the effects are sig- 
nificant only within the widths of the two autoionizing 
states. In this subsection, taking the same test system as 
in Sec, nil Cl we examine the electron and the absorption 



spectra with respect to the intensity and the detuning of 
the laser. In the following, Rabi oscillation is quantita- 
tively express by Rabi frequency 

m = ^J[DEoit)f+S^ (15) 

where D is the dipole matrix element between the cou- 
pled states, Eo{t) is the envelope of the pulse, and S is 
the detuning of the carrier frequency. For near-resonance 
condition, pulse area given by 

/oo 
n{t)dt, (16) 
-oo 

represents the population transfer between the states at 
the end of the pulse; a thorough transfer from one to the 
other is indicated hy A = n. The numerical integration 
covers the whole temporal range of the coupling pulse. 

Figure [5] shows the calculated spectra for a fixed time 
delay to = 15 fs and various laser intensities. For a typi- 
cal two-level system, each cycle of Rabi oscillation revives 
the initial state and changes the phase of its coefficient by 
TT. Similarly, for the autoionizing states here, the bound- 
state coefficients change the phase by it after each cycle. 
For the lowest laser intensity, i.e., Il = 0.5 TW/cm^ and 
A = 1.347r, the inverse Fano line shape (the g-parameter 
changes sign) in Fig. [SJa) and the single peak in (b) 
are formed where a fraction of the wave packet carry- 
ing the changed phase swings back to 2s2p. This con- 
dition has been discussed in details in Ref. |2^. For 
Il = I TW/cm^ or A = IMn, the electrons are ex- 
pressively driven back to 2s2p. The inverse Fano peak 
is prominent in Fig. [Sja) and (c), and the 2s^ profile 
is depleted by about half in Fig. Mjo). As 1^ further 
increases, e.g., A = 2.677r and 3.787r for II ~ 2 and 
4 TW/cm^ respectively, in the same laser duration the 
two levels undergo more oscillatory cycles before autoion- 
ization, which results in more interference fringes within 
the covered energy region. 

The XUV absorption spectra in Fig. ^c) are, as ex- 
plained earlier, almost the same in shape as the 2s2p elec- 
tron profile but with stronger signals. This implies that 
the XUV absorption spectra would yield essentially the 
same information as that from the electron spectra. Since 
high-resolution spectra are harder for electrons than for 
photons, the detailed structures predicted here are bet- 
ter carried out using photoabsorption spectroscopy. On 
the other hand, as shown in Fig. [5jc) , at larger II , the 
absorption signals near 60.3 eV appear to be negative; it 
means some amount of energy in that frequency range is 
released by the atom in the form of XUV emission. In 
Fig. [5ljd) we further note that the laser absorption spec- 
tra feature an oscillatory pattern with the same nodal 
points in energy. When the laser is intensified, the nodal 
points do not shift, but the amplitude of the pattern in- 
creases. This suggests that the absorption and the emis- 
sion of the laser simultaneously increase. 

So far we have discussed the special resonance con- 
dition where the 540 nm laser resonantly couples 2s2p 
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FIG. 5: (Color online) Photoelectron spectra near the (a) 
2s2p and (b) 2s^ resonances, and photoabsorption spectra 
for the (c) XUV and (d) laser fields, for different peak laser 
intensities (labeled in the units of TW/cm^) and fixed time 
delay to = 15 fs. The gray curve is the spectrum measured 
without the laser. The horizontal black line indicates zero 
signal. 

and 2s^. If the wavelength shifts, the detuning should 
be taken into account, and the spectra are expected 
to change accordingly. In Fig. [51 two additional wave- 
lengths are plotted, where the detuning is about 0.8 eV 
above and below the resonance, respectively, with fixed 
II = 0.7 TW/cm^. Due to the detuning, the domi- 
nantly populated electron energy of the 2s^ resonance 
is shifted upward or downward as the laser photon en- 
ergy is lower or higher respectively, thus shifting the 2s^ 
profile as shown in Fig. [SJb). The laser absorption pat- 
terns for all three wavelengths in Fig. Hl^d) are formed 
by the same interference mechanism but different carrier 
frequencies, which consequently show overall shifts in en- 
ergy but the same geometric features. The 2s2p profile 
in Fig. HJa) and the XUV profile in (b) arc obviously de- 
formed by the detuning. However, a strong indication of 
the 2s2p resonance energy, at 35.55 eV in electron en- 
ergy or 60.15 eV in XUV energy, is unaffected. With 
the 0.8 eV detuning, its is the wing of the incident laser 
profile that couples the two states, where the slope of the 
wing pickes up additional phase and strength dependence 
in the Rabi oscillation, which then further deforms the 
spectra from the 540 nm case. 



IV. SUMMARY AND CONCLUSION 

We have developed a theoretical model to calculate the 
photoabsorption spectra of an atom for an attosecond or 
femtosecond XUV pulse in the presence of a moderately 
intense laser pulse. The XUV pulse pumps the atom 
to an autoionizing state when a time-delayed laser pulse 
couples this state to another autoionizing state. The 



FIG. 6: (Color online) As Fig. [5] but for different laser wave- 
lengths (labeled in nanometers) and fixed peak laser intensity 
II = 0.7 TW/cm^ and time delay to = 15 fs. 

model is extended from the standard theory of KIT to 
two autoionizing states and to ultrashort pulses. The 
wavcfunction in this problem was formulated and ap- 
plied to calculate photoelectron spectra previously [26l |. 
In this paper, we generalize the theory to calculate pho- 
toabsorption spectra. To illustrate that the model is 
applicable to both attosecond and femtosecond pulses, 
the time-delayed absorption spectra of two available re- 
cent exp eriments are simulated, one using SAPs on ar- 
gon |27l | and the other using APTs on helium [2l[. For 
both cases, good general agreement with experimental 
data have been found. However, in both experiments, the 
coupling IR laser severely ionizes the autoionizing states, 
and the influence of laser coupling on the resonances is 
less spectacular. Furthermore, the spectral resolution in 
both experiments are still limited and unable to reveal 
the rich structures of the IR-modified autoionizing pro- 
files. 

To highlight the effect of laser coupling, we suggest to 
use a 540 nm laser to resonantly couple the autoioning 
2s2p and 2s^ states in helium. Because of the large bind- 
ing energies, ionization by the dressing laser is negligible 
for both states, and the coupling effect is maximized. We 
then investigate the electron profiles of both resonances 
and the absorption profiles of both the XUV and the 
laser pulses. The result shows that the XUV absorption 
profile near 2s2p is almost the same as the correspond- 
ing electron emission from the 2s2p resonance, but the 
absorption signal is higher than the electron signal. The 
absorption and electron profiles within the widths of the 
resonances reveal the resonance features modulated by 
the dressing pulse. For such studies, photoabsorption 
measurements are more favorable since better spectral 
resolution can be reached with photons. 

We should emphasize that the present calculations 
have been limited to light interactions with single atoms 
only. For photoabsorption spectroscopy, light pulses are 
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directed in the gas medium coherently and thus the effect 
of macroscopic propagation should be included [2^ . This 
would offer a practical method to shape the XUV pulse 
by supplying an time-delayed intense IR laser in a gas 
medium, where the propagation effect can be modeled 
after the high-order harmonic generation in gases [s^ . 
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